Irrelevance of information outflow in opinion dynamics models 
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The Sznajd model for opinion dynamics has attracted a large interest as a simple realization of 
the psychological principle of social validation. As its most salient feature, it has been claimed that 
the Sznajd model is qualitatively different from other ordering processes, because it is the only one 
featuring outflow of information as opposed to inflow. We show that this claim is unfounded by 
presenting a generalized zero-temperature Glauber-type of dynamics which yields results indistin- 
guishable from those of the Sznajd model. In one-dimension we also derive an exact expression for 
the exit probability of the Sznajd model, that turns out to coincide with the result of an analytical 
approach based on the Kirkwood approximation. This observation raises interesting questions about 
the applicability and limitations of this approach. 
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In the last decades statistical physics has crossed many 
boundaries between different fields, becoming, with its 
methods and concepts, a powerful tool for the inves- 
tigation of a broad range of disciplines. This process 
has been mutually beneficial, since the consideration of 
problems far away from a purely physical motivation 
has greatly broadened the kind of theoretical questions 
and conceptual challenges statistical physics is called to 
tackle. One of the settings in which this cross-fertilization 
has been particularly fruitful is opinion dynamics 
where the goal is to understand how global consen- 
sus/understanding/agreement emerges out of disorder, 
based on local interactions. In this field many simple 
models akin to those of statistical physics have been in- 
troduced, both by social scientists and physicists 
leading to intense activity and remarkable results. In this 
context, the model introduced by Sznajd- Weron and Sz- 
najd 0, commonly denoted as Sznajd model (SM), has 
enjoyed an exceptional success as the first one encod- 
ing the principle of "social validation", stating that the 
convincing power of an individual is greatly enhanced if 
another individual supports the same view. 

The dynamics of the SM in one dimension is defined 
as follows [1, 0] : Each site in a one dimensional lattice 
is endowed with a binary variable (spin) = ±1. At 
each time step a pair of neighboring sites is selected at 
random, i and i + 1. If these individuals have the same 
opinion, ai — di+i = cr, the opinion of all the neighbors 
of i and i 4- 1 changes to the common value ct; other- 



wise, nothing happens l21|. The process is iterated until, 
on a finite system, a final consensus (all spins equal) is 
reached. Generalizations to higher dimensions have been 
introduced and are described below. Typical quantities 
of interest are the consensus (fixation) time T(x, N), de- 
fined as the time needed to reach the state with all spins 
equal for a system of size N, starting from a configuration 
with a fraction x of positive spins, and the exit probabil- 



ity E{x), defined as the probability that the final state 
will be all cr,; = +1. 

The Sznajd model is similar to other simple models for 
dynamics of Ising spins in the absence of bulk noise, such 
as the voter model and the zero temperature Glauber 
dynamics. However, much emphasis has been put [13- 



12| on the claim that SM is fundamentally different be- 
cause it is the only model where "information flows out" 
(i.e. spins propagate their state to their neighbors) as op- 
posed to other models where a central spin adapts itself to 
the state of the surrounding ones ("information inflow" ). 
This claim is mainly supported by the shape of E(x) in 
one dimension, which is linear for the Glauber zero tem- 
perature dynamics (as well as for the voter model) , while 



it is nontrivial for SM [13|, |lj| . Also the consensus time 
T{x, N) has a dependence on x for SM that is not found 
in other types of dynamics. 

In this Letter we show this claim to be unfounded, by 
presenting two clearly "outflow" and "inflow" dynamics, 
given by simple extensions of the SM and the Glauber 
models, respectively, in which the number of sites in- 
volved in a single spin update is a model parameter. The 
analysis of these models allows us to check that the pos- 
tulated difference between "inflow" and "outflow" dy- 
namics in fact docs not exist. In particular, we show 
that in one dimension the exit probabilities and consen- 
sus times of both models are the same. The consider- 
ation of the two dimensional and mean-field cases adds 
additional strength to our result. Additionally, we pro- 
vide an exact expression for the exit probability of SM in 
one dimension, revealing that previous results based on 
a Kirkwood approximation are also exact, due to some 
surprising cancellation of errors that remains to be un- 
derstood. 

The models we consider arc defined in one dimension 
as follows: 

Sznajd Model of range R, SM{R): At each time step, a 
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-- Eq.(l) 
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FIG. 1: Main: Exit probability for the SM{R) on a one di- 
mensional lattice of size L — 1000. The full line corresponds 
to the analytical prediction, Eq. ([T}. The number of real- 
izations of the dynamics is M = 10^. Inset: Difference be- 
tween the numerical results for SM{1) and the theoretical 
as L changes. The uncertainty due to sampling error for 
L = 1000 is given by ±^{E{x)[l - E{x)]/M), where the 
number of realizations of the stochastic process is M = 10^. 



pair of nearest neighbor sites, i, i + l, is chosen at random. 
If they share the same state, ai = ct^+i = a, then the 
2R neighbors, to left and right, respectively, change their 
value to CT, i.e. a-j — > a, for j G [i—R,i — l][j[i+2,i+l+R]. 
Otherwise, nothing happens. In this "outflovir" dynamics, 
the opinion of two adjacent equal spins thus extends to 
all their 2R neighbors, the case R = 1 corresponding to 
the standard SM. 

Zero temperature Glauber dynamics of range R, G{R): 
The elementary step consists in randomly selecting a site 
i and evaluating the local field given by the sum of the 2R 
spins in the interval [i~R,i — l]\j[i + l,i + l + R]. If the 
local field is positive or negative, the variable ai aligns 
with it. Otherwise the spin is randomly set to ±1 with 
probability 1/2. For any R the dynamics is obviously of 
"inflow" type, as the central spin is affected by the state 
of surrounding spins. The case G{1) coincides with the 
usual zero-temperature Glauber dynamics. 

Let us consider uncorrelated initial conditions in which 
each vertex has a probability x to be in the -1-1 state, and 
correspondingly a probability 1 — a; to be in state —1. As 
in other ordering processes of this kind, the evolution in 
both SM{R) and G{R) models proceeds in two separate 
stages. Initially, homogeneous domains of spins up or 
down quickly form at small scale. This stage lasts for a 
time of the order of a few Monte Carlo steps per spin. 
Later on, domain boundaries diffuse around and annihi- 
late upon encounter, leading to larger and larger domains 
and eventually to consensus. The duration of this second 
stage grows with the system size as L^. While in the 
first stage the dynamics depends on the model's micro- 
scopic details and the magnetization is not conserved, the 
second regime is very similar for voter, generalized SM 



or generalized Glauber models, with marginal variations 
due only to the details of the annihilation process. In 
this stage, the diffusion-annihilation boundary dynamics 
leads to the conservation of the average magnetization. 

We first study the behavior of the SM{R) model, plot- 
ting in Fig. [1] the exit probability for this model, com- 
puted numerically for different values of R. Remark- 
ably, E{x) turns out to be completely independent of 
the range of the interaction R. By taking advantage of 
this independence of R, we can derive the exact form of 
the exit probability, which is very easy to compute for 
R> [L — 2)/2. In such a case, the diffusive regime is ab- 
sent and the system becomes fully ordered after the first 
successful microscopic update. The dynamics proceeds 
by choosing at random 2 consecutive sites, that will be 
both in state +1 with probability x^, in state —1 with 
probability (1 — x')^, and in a mixed state with probabil- 
ity 1 — a;^ — (1 — a;) ^. The exit probability is given by 
the probability that a pair of sites in state -|-1 is chosen 
before any pair of sites in state —1. So, we can write 



E{x) = a;2 ^[1 - a;2 - (1 - a;)^]" = 



n=0 



a;^ + (1 — xY 



(1) 



Another way to derive Eq. ([T]), valid for smaller values 
of R, is as follows: In the initial stage each successful 
update will give rise to a domain of 2 -I- 2R equal sites, so 
that in a time of order unity the system will be roughly 
subdivided into L/{2+2R) domains of size of order 2-|-2i?. 
At the end of this stage [l^l, the density of +1 spins 
wiU be x' = 1 X x^/[x^ -f (1 - a:)2] + x (1 ~ x)^/[x^ + 
(1 — a;)^] . In the ensuing second stage the conservation of 
magnetization implies that the exit probability is E{x') = 
x\ independent of domain size, yielding again Eq. ([1}. 

Fig. [1] shows that Eq. ^ provides a very accurate de- 
scription of the exit probability of the generalized SM. 
The inset of the figure proves moreover that the small 
deviations of the numerical results for E{x) around the 
theoretical value can be fully ascribed to fluctuations 
around the expected value due to the finite number of 
realizations of the process. This confirms that Eq. ([T]) is 
the exact solution of the exit probability for the SM. 

It is crucial to remark that Eq. ([T|) coincides with the 
expression for the exit probability of SM calculated by 
solving analytically the hierarchy of equations for multi- 
spin correlation functions within a Kirkwood approxi- 
mation decoupling scheme [l^, [l3|- In this case then, 
the Kirkwood approximation turns out to provide an ex- 
act solution for the SM model. This is indeed a striking 
result, since numerical tests show that the assumptions 
made in the Kirkwood approximation are largely violated 
during the dynamics. 

Turning now to the generalized Glauber dynamics 
G{R), numerical simulations for R ~ 2 (Fig. [2]) show 
that also in this case E{x) is in excellent agreement with 
Eq. ([Ij. Hence the exit probability of SAI{R) and of the 
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FIG. 2: Exit probability for the G(2) dynamics in one dimen- 
sional lattices of increasing size. The full line corresponds 
to the analytical prediction for the SM, Eq. ([T]). The num- 
ber of realizations of the dynamics is 10®. Inset: Numerical 
values of the consensus time T{x) rescaled by its maximum 
value T{x = 0.5), for G(2) and SM{R), compared with the 
analytical prediction (see text). Lattice size L = 500. 



G(2) model arc indistinguishable. The closeness of the 
two models is further confirmed by the inset of Fig. [21 
where the consensus time T{x) (divided by T{x = 0.5) 
to factor out trivial temporal rescalings) is reported: The 
time needed to reach the final consensus is the same for 
both SM{R) and G{R = 2) models. Fig. [2] provides fur- 
ther evidence of the independence of the generalized SM 
with R and allows to conclude that direction of "infor- 
mation flow" is irrelevant: The behavior of Sznajd model 
with "outflow" dynamics coincides with the behavior of 
G{R = 2) model, based on "inflow" . 

The dynamical division in two stages, illustrated above 
to derive the exit probability, is useful also to obtain 
an analytical estimate of the time T{x) to reach con- 
sensus for the SM. As described above, in a time of order 
unity the density of -1-1 spins reaches its asymptotic value 
x' = x^ l\x^ + (1 — a;)^], where x is the magnetization in 
the initial state. The ensuing evolution is essentially the 
same followed by the voter model, for which the consen- 
sus time is known and whose dependence on the initial 
density of up spins x' is T oc —\x' lna;'-|-(l — x') ln(l — x')]. 
Expressing x' in terms of the initial value x one obtains 
an analytical formula for the consensus time of the SM. 
The comparison with numerics (Fig. [21 inset) is rather 
good, the discrepancy observed being probably ascrib- 
able to the slightly different behavior of the models when 
two boundaries are one site far apart. While in the voter 
dynamics they have equal probability to collide or to go 
to distance 2, they deterministically collide in SM. 

The strong relationship between the 0(2) and SM is 
not limited to one-dimensional systems. Let us consider 
a random neighbor topology, i.e. a fully connected sys- 
tem where the interaction occurs with neighbors chosen 
randomly at each time step. Slanina and Lavicka [14 1 



FIG. 3: Consensus time T(a;, A'^) as a function of the initial 
fraction x of -f 1 spins, for SMiX) and G(2) models on a ran- 
dom neighbor topology. Upper curves are for size A*" — 10*. 
Lower curves are for size A'^ = 10'^. 



have analyzed the standard Sznajd dynamics in this case, 
characterized by the transition rates 

x^[\-x), 



Prob[a; 



1/A^] 



Prob[a; a; - l/AT] = x{\ - xf 



(2a) 
(2b) 



where N is the system size. For the generalized Glauber 
G{2) dynamics the rates can be also easily worked out 

Prob[a: -> x + 1/iV] = x^{\- x)(i-2x), (3a) 
Prob[x a; - 1/iV] = x(l - a;)2(l + 2x). (3b) 

The only variation is given by correcting factors which 
are smooth and positive, thus implying that no basic fea- 
ture of the dynamics will change. In particular, following 
the inverse Fokker-Planck formalism [Hj, it is possible 
to show that the exit probability takes in both cases the 
form of a Heaviside step function, E{x) = 0(a; — 0.5) for 

— > cxD, as expected due to the presence of an imbalance 
between the rates in Eqs. ([2]) and Concerning the 
consensus time, we report the results of computer sim- 
ulations in Fig. [21 which prove the equivalence between 
the SM and the G(2) models at the mean-field level. 

In finite dimensions larger than one, many possible 
ways to define SM have been introduced Simi- 
larly, there are various possibilities to define the Glauber 
dynamics for generic R. We select the following ones. 
For the G(i?) model, the local field for a site (i,j) is 
given by the sum of all spins up to the i?-nearest neigh- 
bors. In particular for i? = 2 in c? = 2 the local field is 
given by the sum of the eight spins surrounding (i, f) and 
forming together a square of side 3. For the SM, on the 
other hand, we consider two variants. In SM-I dynamics 
a bond is randomly chosen, either along the vertical or 
horizontal direction, and if the sites at the extremes of 
the bond are equal, all the 6 nearest neighbors of both 
sites are updated accordingly. In SM-II, we select a pla- 
quette of four sites and if they are in same state, the 8 
nearest neighbors are made equal. 
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FIG. 4: Density p{t) in d = 2 for SM-I, SM-II, G{1), and G(2) 
dynamics. System size L = 10000. Inset: p(t) for several runs 
on a SM-I system of size L = 1000, illustrating the stripe 
phenomenon. 



The probabiUty E(x) to end up with all +1 spins is 
for all variants of SM given (in the large size limit) by a 
step function E{x) = e{x - 0.5) ^. As can be expected 
based on the fact that the dynamics is driven by curva- 
ture 



17[, the same occurs for the G(2) modeL provided 
no freezing in a striped configuration occurs jl8| . This 
phenomenon, which affects asymptotically G(l) dynam- 
ics in (f = 2 with a finite probability [18[ (and clearly 
affects G(2) dynamics as well), is present also in the evo- 
lution of the SM. In this case, straight stripes along one 
direction in a two-dimensional lattice are not fully stable, 
given the intrinsic destabilizing mechanism present in Sz- 
najd dynamics at microscopic scales. Nevertheless stripes 
do often form during the evolution and they persist for 
very long times. The presence of stable or long lived 
metastable striped states makes a comparison between 
the consensus time T{x, N) in SM and Glauber G{R) 
models impossible. A quantity allowing a better analy- 
sis of the ordering behavior of two-dimensional systems 
is the fraction p{t) of nearest neighbor pairs that are in 
opposite states. Figure |4] shows that for Sznajd and gen- 
eralized Glauber dynamics (with both R = 1 or R = 2) 
the evolution is the same, apart from irrelevant transients 
and global temporal scales: The density p decreases as 
the signature of curvature-driven coarsening dy- 
namics [l3|- On the other hand, the plateaus exhibited 
in some realizations of Sznajd dynamics for long times 
(Fig.m inset) indicate the effective presence of long-lived 
metastable states. The perfect analogy between Sznajd 
and Glauber dynamics goes beyond the decay of p{t). 
The scaling functions for the two-point correlation func- 
tion C{r,t) (not shown) are virtually the same. 

In summary, we have shown that the behavior of Sz- 
najd model for opinion dynamics has no feature that 
distinguishes it from a generalized zero temperature 
Glauber dynamics for Ising spins. While in dimension 
d > 1 this could be expected on the basis of general 



considerations on coarsening systems, in d = 1 this re- 
sult is highly nontrivial. In one-dimensional systems, the 
standard Sznajd dynamics actually differs from the usual 
zero-temperature Glauber dynamics, as it has been ex- 
tensively reported in the literature. However, when the 
range of the interactions is extended to i? = 2, the gen- 
eralized Glauber dynamics is indistinguishable from Sz- 
najd. The conclusion is that "outflow" dynamics is not 
qualitatively different from "inflow" dynamics. A pos- 
sible objection to this conclusion is that inflow and out- 
flow dynamics are actually different because SM{R) docs 
not depend on R, while G{R) dynamics does. This ar- 
gument is rebutted by considering another extension of 
SM, in which the number of equal spins needed to con- 
vince neighbors is a parameter q. Numerical and ana- 
lytical arguments, to be reported elsewhere [l^, show 
that q strongly affects the dynamics and that such a gen- 
eralized "outflow" dynamics gives results very close to 
those of the G{R) "inflow" model with R = q. While 
studying the equivalence of SM{R) and G(2) we have 
derived an exact formula for the exit probability of SM, 
that turns out to coincide with the one obtained by using 
a Kirkwood approximation. Notice that, on the contrary, 
Kirkwood approximation fails for the G(2) dynamics [19| . 
These findings call for additional research to understand 
when the Kirkwood approximation works, when it fails, 
and how it can be systematically improved. 
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